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Recent studies of strongly interacting atoms and photons in optical cavities have rekindled interest in the 
Dicke model of atomic qubits coupled to discrete photon cavity modes. We study the multimode Dicke model 
with variable atom-photon couplings. We argue that a quantum spin glass phase can appear, with a random 
linear combination of the cavity modes superradiant. We compute atomic and photon spectral response functions 
across this quantum phase transition, both of which should be accessible in experiments. 
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Introduction. Ultracold atoms in optical cavities have 
emerged as attractive new systems for studying strongly- 
interacting quantum many body systems. Photon exchange 
can mediate long-range interactions between the atomic de- 
grees of freedom, and this opens up rich possibilities for cor- 
related phases. In the celebrated atomic realizations of the 
superfluid-insulator quantum phase transition [1|, the light 
field acts in a secular manner, creating a potential which traps 
the atoms in an optical lattice; consequently the atom-atom 
interactions are only on-site, and this limits the range of pos- 
sible phases. In contrast, the seminal recent experiments of 
Baumann et al. [2,3], realizing a supersolid phase, have long- 
range interactions mediated by active photon exchange J4] . 

Baumann et al. argued that their experiments could be de- 
scribed by the Dicke model, as in the proposal of Nagy et 
al. 0. The Dicke model couples photons in a single cav- 
ity mode uniformly to N atomic two-level systems ('qubits'). 
In the limit — > oo, this model exhibits a phase transition 
||6rl9l to a "superradiant" phase when the atom-photon cou- 
pling is strong enough. In terms of the qubits, the superra- 
diant phase is a 'ferromagnet' which spontaneously breaks a 
global Ising symmetry, and so we refer to it as FMsr. In the 
experiments by Baumann et al, the superradiance of the cav- 
ity photon mode is accompanied by 'self-organization' of the 
atoms into a density wave pattern |[T0ffl"2l . 

Here we study extensions of the Dicke model to multiple 
photon cavity modes, and with non-uniform couplings be- 
tween the atomic qubits and the photon modes. Spatial mode 
variations for the single-mode Dicke model were considered 
in Ref. 13 Multimode Dicke models have been studied earlier 
(6] [14l4T6ll . but were simplified by ignoring the variations in 
the atom-photon couplings. We argue here that qualitatively 
new physics emerges in the multimode case when the spatial 
variation is treated seriously. We show that large variations 
in the atom-photon couplings can give rise to a quantum spin- 
glass (QSG) phase. We will describe quantum-critical dynam- 
ics associated with the onset of this spin glass order. 

Dimer et al. IfTTl have discussed an experimental realiza- 
tion of the Dicke model using internal atomic degrees of free- 
dom, that is, Raman transitions between multiple atomic lev- 
els. We expect that such schemes can be generalized to a mul- 
timode Dicke model that respects a global Ising symmetry, 
which is then spontaneously broken in the FMsr and QSG 



phase, respectively. More specific realizations of the multi- 
mode Dicke model were described recently by Gopalakrish- 
nan et al, in a paper [18 1 which appeared while our work 
was being completed. The same authors had previously out- 
lined how Bose-Einstein condensates in multimode cavities 
can lead to frustration and glassy behavior IfTTl IT21 . Such 
experiments on the multimode Dicke model would provide 
a unique realization of a quantum spin glass with long-range 
couplings, and provide a long-awaited testing ground for theo- 
ries of quantum systems with strong interactions and disorder. 
Condensed matter realizations of quantum spin glasses have 
shorter-range couplings, and so do not directly map onto the 
theoretically solvable systems analyzed in the present paper. 

Before describing our computations, we point out a key 
distinction between the transitions involving onset of FMsr 
versus QSG order. In the single-mode Dicke model, all the 
qubits align in a common direction near the FMsr phase, and 
can therefore be described by a collective spin of length N/2 
which behaves classically in the limit of large N. Conse- 
quently, the dynamics near the phase transition can be de- 
scribed by classical equations of motion [19], and the single- 
mode Dicke model does not realize a quantum phase transi- 
tion in the conventional sense of condensed matter physics. In 
contrast, we will argue here that the onset of QSG order in the 
multimode Dicke model has non-trivial quantum fluctuations 
even in the limit of large N, and the critical properties cannot 
be described by an effective classical model. Experimental 
studies are therefore of great interest. 

Model. The Hamiltonian of the multimode Dicke model is 

M . N N M 

H = Yj w < fl i »<■ + T Yj °"c + Yj Yj 8it ( ai + fl <' ) a t- (l) 
i=i e=\ e=i ;=i 

This describes N two-level atomic qubits with level splitting 
A/2 and M photon modes with frequencies w, coupled by an 
atom-photon coupling g% which depends on the photon (z) and 
atom ({) number. a\, a,- are bosonic creation and annihila- 
tion operators, respectively, fulfilling canonical commutation 
relations. cr^ z are spin- 1/2 operators with Pauli matrix rep- 
resentation. As explained in detail in Refs. IfTTl [T8l . the two 
states of the Ising spin in Eq. (JTJ) map onto two different stable 
ground-state sublevels, |1) and |0), of three-level A atoms. |1) 
and |0) are indirectly coupled through a pair of Raman transi- 
tions to an excited state \e) which are driven by the classical 
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field of a pair of external lasers. Upon adiabatic elimination 
of \e), one obtains Eq. (JTJi with <x~ = |lf)(lf| - |0f)(0^| and 
°~( - |lf)(0f| + |0f)(lf|. The parameters <w;, A, and gn can be 
controlled through laser frequencies and intensities. This tun- 
ability enables access to the strong-coupling Dicke regime. A 
dispersive shift of the cavity frequencies ~ a\aj<T z does not 
modify our results significantly, and so will be set to zero 
for simplicity. A simple choice for a spatially varying atom- 
photon coupling is g^ = g cos (kjX[) with k, the wavevector of 
the photon mode, and x c the coordinate of atom t. 

In the single-mode, large photon wavelength case, we have 
M = 1, u>i = wo, and gn = gj V^V and the model can be 
solved exactly in the N — > oo limit EG). At zero temperature, 
there is a continuous phase transition between a paramagnetic 
phase (PM) and a superradiant ferromagnetic phase (FMsr) 
at g — g c = VAwo/8 at which the Ising symmetry (a, cr x ) — > 
(-a, -<x v ), is spontaneously broken. 

For the multimode Dicke model, it is useful to integrate out 
the photon degrees of freedom in a path-integral representa- 
tion. Then the qubits are described by a Hamiltonian similar 
to the Ising model in a transverse field, 

A N 1 

£=1 tm 

The exchange interactions 7f m are mediated by the photons 
and have a frequency dependence associated with the photon 
frequencies w,; thus Eq. (j2j) is to be understood as an action 
appearing in an imaginary time path-integral summing over 
time-histories of the qubits. The long-range exchanges 

M - 

jum = Z . (3) 

depend on Q., the imaginary frequency of the qubits in the 
path integral. Note that although we have formally integrated 
out the photons, we demonstrate below that the photon-photon 
correlation function is directly related to the atom-atom cor- 
relation function as obtained by solving Eq. (j2j). 

If we ignore the frequency dependence in Eq. (BJ, the J( m 
have a structure similar to the Hopfield model of associative 
memory 11201 , with M 'patterns' g^. For M small, it is ex- 
pected that such a model can have M possible superradiant 
ground states with FMsr order (o~f} K git, i = 1 • • • M. In 
the spin-glass literature, these are the Mattis states which are 
"good" memories of the patterns g [20 1 . The critical proper- 
ties of the onset of any of these FMsr states should be similar 
to those of the single mode Dicke model. 

Our interest in the present paper is focussed on larger values 
of M, where the summation in Eq. ([3]l can be viewed as a sum 
over M random numbers. Then, by the central limit theorem, 
the distribution of Jf,„(Cl) is Gaussian. Alternatively, the ran- 
domness of J( m (Gl) can be enhanced by passing the trapping 
laser beams through diffusers so that atomic positions X( are 
randomly distributed inside the cavity Ifl8l . In either case, we 
assume a Gaussian distribution characterized by its mean and 
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FIG. 1: (Color online) Zero-temperature phase diagram for ojq = 1, 
A = 1 computed from Eq. yTj . PM means paramagnet, FM S r super- 
radiant ferromagnet, and QSG quantum spin glass, q is the Edwards- 
Anderson order parameter and if/ is the atomic population inversion 
or ferromagnetic order parameter. 

variance 



JeJm = MO)/N (4) 

6 J J = (See 6 em .)K(n,Q.')/N, 

where the line represents a disorder average, and 6J( m is the 
variation from the mean value. We have assumed couplings 
between different sites are uncorrelated, and this will allow 
an exact solution in the — > oo limit, modulo an innocuous 
softening of the fixed length constraint on the Ising variable 
ET1I221 . We will allow arbitrary frequency dependencies in 
Jq(Q) and K(Q., £2'). The factors of N ensure an interesting 
N — > 00 limit [23 1. Especially for finite M, one could also use 
the methods of Ref. [20 1 to extend our analysis to models in 
which the gn rather than the J( m {£X) are taken as independent 
random variables. However, as long as the photon modes can 
be chosen so that the 7f m (Q) vary in sign and magnitude, our 
analysis should remain qualitatively correct also for smaller 
values of M. 

Key results. We will show below that, in the limit of large 
atom number N, the results depend only upon 7o(^ = 0) and 
K(Cl, -Q). Here, we will display the phase diagram and spec- 
tral response functions for the simple choices Jo(0) = 2g 2 /a>o 
and K(Q, -Q) = 7 2 (Q) with 

J(D.) = 2f 2 wo/(Q 2 + Wq). (5) 

In Fig. [TJ we depict the ground state phase diagram; a re- 
lated phase diagram in a condensed matter context was ob- 
tained in Ref. |24] All phase transitions are continuous and 
the respective phase boundaries merge in a bicritical point at 
(& = 0-086,4 = 4). 

The intersection of the PM-FMsr phase boundary with the 
vertical axis at t 2 = corresponds to the phase transition in 
the single mode Dicke model without disorder |[8][T7). In this 
case, a number of analytical results can be obtained from Eq. 



3 



FM 



SR 



0.0 0.2 



co = 1.3 
w = 1.015 

to = 0.5 



n 



ca 
a 

V 2 



i QSG 




O) = 

U/Q — 


1.7 
1.59 


// \\ 

/ \ 

/' \ 
/' \ 

/' V 
/' * 
II 


/ 

\ / 


CO Q = 


0.7 


h 
a 

li i 


i\J 

T\ 
\ \ 






j \. 


\ \ 

\ \ 
\ \ 
1 x 




1 



1.0 

n 



FIG. 2: (Color online) rf spectral response function of the atomic 
qubits in the FMsr phase for various photon frequencies and t 2 = 
0.025, g 2 = 0.2, A = 1. Q is a real measurement frequency. The 



red arrow at Q, 
weight q ~ if/ 2 from Eqs 

above Eq. ql3p. For the Dicke model without disorder (/ 



illustrates the delta function contribution with 
13} . The value of the gap is given 
0), the 
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spectral function following from Eq. |6]( consists of nothing but two 
delta functions: the red arrow at Q = and the grey arrow at CI = 
g V2A/o> (plotted for cj = 1.015). 



(Hi, in agreement with the earlier work. The critical atom- 
photon coupling is g 2 - Awo/8 and the local cr x ( spin suscep- 
tibility in the FMsr phase is (for imaginary frequencies) 



er(")U 



Q 2 + 2Ag 2 /to 



+ (A 2 2tt6(Q.). 



(6) 



The corresponding radiofrequency (rf) spectral response 
function of the atomic qubits for real frequencies, 
-Im[g flfl (/Q -> Q + /0+)], is depicted in Fig. [2] The su- 
perradiance is encoded in the zero frequency delta function 
contribution, whose weight is proportional to the atomic pop- 
ulation inversion if/. However, away from the zero frequency 
delta function, there is a spectral gap, and the remaining 
spectral weight is a delta function at frequency ^2Ag 2 /uj a . 

The superradiance also appears as a photon condensate 
(a,) = - 'Z,e(gic/(2cL>i))(cr i e ;). We have computed the atomic 
population inversion, {erf) = if/, and the Edwards-Anderson 



order parameter (cr x ( ) 2 = qqsg in Eqs. ( 13 



14 1. Both of these 



are related to (a,), but computation of the latter requires more 
specific knowledge of the gu. For Q + 0, the photon correla- 
tion function follows from Eq. ([6]) 



(iQ-cjdSij + YjgicgjcQTi^) 



(7) 



where the right-hand-side is a matrix inverse, as can be ob- 
tained from integrating out the atomic fields from the path- 
integral representation of Eq. ([TJ. 

Upon introducing small disorder (with t + 0), as long as we 
remain in the FMsr phase, the zero frequency delta function 
and spectral gap survive, although the higher frequency spec- 
tral weight changes, as shown in Fig. [2] This spectral gap is 



FIG. 3: (Color online) rf spectral response function of the atomic 
qubits in the QSG phase for various photon frequencies and t 2 = 
0.175, g 2 = 0.05, A = 1. The red arrow at Q = illustrates the delta 
function contribution with weight q ~ <5>q SG from Eqs. l |12|14[ l. 



present across the phase transition from the FMsr phase to the 
PM phase. Thus all the low energy fluctuations in the critical 
theory for this transition are restricted to the zero frequency 
delta function, which can be described in classical theory for 
the spins: this is the reason this transition is more properly 
considered as a classical phase transition. 

For a sufficiently large value of f 2 , the system undergoes a 
quantum phase transition to the QSG ground state. In con- 
trast to the PM-FMsr transition, at the QSG transition, and in 
the entire QSG phase, there is spectral weight at a continuum 
of frequencies reaching zero (see Fig. |3J. Thus the onset of 
QSG order from the PM phase is a genuine quantum phase 
transition, whose universality class was described in Refl22l 

The PM phase is clearly delineated from both, the QSG 
and the FMsr phases: the PM phase has a gapped spectral 
response and no superradiant photon condensates. 

We also note that in all phases, while the spectral function 
has a universal form at low frequencies, its high frequency 
behavior is strongly dependent upon the forms of Jq(Q) and 
J(£l). For the forms in Eq. ([5j, the spectral function is sup- 
pressed to zero at Q = ojq. 

Experimental signatures. The rf spectral response function 
of the atomic qubits presented in Figs. 2[3 should be observ- 
able via radiofrequency spectroscopy 1251 1261 . 

Measuring the spectrum of photons leaving the cavity 
through its imperfect mirrors at loss rate k allows for an in- 
situ measurement of our phase diagram, Fig. [T] Our predic- 
tion for the spectrum of intra-cavity photons, Eq. d7), can be 
related to the extra-cavity photons via the input-output formal- 
ism lfT7l l27l l28l . For this case of a dissipative Dicke model, 
we note a similarity of the decay effects to those in theories 
of metallic spin glasses |29|, in which the spin qubits are 
coupled to a "reservoir" of continuum spin excitations near 
the Fermi surface. This coupling leads to a damping term in 
the dynamics of each spin, but does not significantly modify 
the spin-spin interactions responsible for the spin glass phase. 
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Similarly, for the dissipative Dicke model, decay into photons 
outside the cavity will introduce various damping terms e.g. 
a k|Q| term in the denominator of Eq. pj. As in the previous 
analyses ||29ll , we expect that the quantum spin glass transi- 
tion will survive in the presence of damping, although there 
will be some changes to the critical properties [ 30 1 . 

As in other glasses, we expect slow relaxational dynamics, 
along with memory and aging effects in the QSG phase which 
should be observable via local spin addressing protocols and 
measuring the spin relaxation time scale lfT8ll . 

Conclusion. Observations of these effects in quantum op- 
tic systems would be remarkable. Moreover, the spin glass 
physics is driven by long-range interactions, and this makes 
the theoretical models analytically tractable. We therefore 
have prospects for a quantitative confrontation between the- 
ory and experiment in a glassy regime, something which has 
eluded other experimental realizations of spin glasses. 

Details of the calculation. As discussed in Refs. |2T] [221 
each Ising qubit, with on-site gap A/2, is conveniently rep- 
resented by fluctuations of a non-linear oscillator </>((r) (r is 
imaginary time) which obeys a unit-length constraint. Their 
action at temperature T is then 



I N ri/T 



(8) 



where t is imaginary time, and the Ac are Lagrange multipliers 
which impose the constraints. The only approximation of this 
paper is to replace the Ac by their saddle-point value, iA( = 
A, and to ignore their fluctuations. For decoupled oscillators, 
this saddle-point value is A - A 2 /4, the <p susceptibility is 
A/(£2 2 + A 2 /4), and the resulting gap, A/2, has been matched 
to that of the Ising spin. 

The interactions between the qubits are accounted for as 
before [21 J: we introduce replicas a — 1 . . . n, average over the 
Jc, n using Eq. Q, decouple the resulting two-</> coupling by 
Hubbard-Stratonovich transformation using a ferromagnetic 
order parameter ^"(Q), and the four-</> coupling by the bilocal 
field e ai (Qi, Q 2 ) G2 (the Q are Matsubara frequencies). The 
complete action is 



S = J]s [<p a ,A a ] + 7-^7 (O) 



aft 

r- 



N 



|*T a (Q)| 



-G a6 (-o,-n')X^(^(^') 



|ie afc (af2')i 2 



(9) 



Now we perform the Gaussian integral over the <pf. the re- 
sulting action has a prefactor of N, and so can be replaced by 
its saddle-point value. By time-translational invariance, the 
saddle-point values of the fields can only have the following 
frequency dependence 

>T(Q) = feo/r)^ 
e fl/ '(Q,Q') = (6 n+ a',o/T)[x(£l)d ab + (6 a ,o/T)q], (10) 



and we take A a = A. We have assumed replica symmetry for 
the Edwards-Anderson order parameter q because our inter- 
est will be limited here to T — where there is no replica 
symmetry breaking Il22l . Now the values of the ferromagnetic 
moment ift, the spin susceptibility ,^(£2), q, and A have to be 
determined by optimizing the free energy. The latter is ob- 
tained by inserting Eq. ( 10 1 in Eq. d9]l; after taking the replica 



limit n — > 0, we have the free energy per site 



7 (0)(A 2 T 

n 



-"M")l 2 



+ -K(0,0) X (0)q 



z-( 



K(0, 0)q + Jl(0)ip 2 



A/A-K(0,0)x{0) 



(11) 



Note that this free energy depends only upon 7o(0) and 
K(£l, -Q), as claimed earlier. Our results described in Eq. |6]) 
and Figs. T[3 are derived from a set of coupled saddle-point 
equations obtained from varying Eq. (Hi with respect to^(O), 
q, and A for every Q.. Subsequently we let T — > 0. 

For the choices for K{Q., -Q) and 7o(0) of Eq. Q, the rf 
spectral response function of the atomic qubits plotted in fig- 
ures |2l3]is given by the expression: 



(12) 



-Im[e flfl (;Q ->O + z'0 + )] = 

ol - Q 2 | ^16A 2 t 4 col -(A- Q 2 ) 2 [u>l - Q 2 ) 2 
8Af 4 w 2 



+ q27T6(Q.) . 



The first term is non-zero only for frequencies D. so that the 
expression underneath the square-root is positive. The value 
of the Lagrange multiplier in the FMsr is pinned to Apm — 
A (7 (0) + ^(0, 0)//q(0)). The value of the gap in Fig. || is 

|^fm + <^o - -^16Af 2 <Do + (/Ifm _ j ■ This expression 

equates to zero in the gapless QSG phase shown in Fig. [3] 
where /Iqsg = 2A V^(0, 0). This gap vanishes logarithmically 
faster than (f 2 -f 2 ) when approaching the QSG phase boundary 
due to the square-root behavior of the spectral weight ||2T|[3T| . 
The ferromagnetic moment obtains as 



J 2 Q (0)-K(0,0) 
J 2 (0) 



_~, Z7T IA=A F mI 



(13) 



and if/ vanishes continuously at the FMsr-QSG phase bound- 
ary (at which 7o(0) = y/K(0, 0)) with exponent jSfm = 0.5. 
As per the discussion above Eq. the corresponding pho- 
ton condensate (a,) vanishes with the same exponent. Note 
that the Edwards-Anderson order parameter q is continuous 
across this transition and in the QSG phase given by: 



qQSG 



1 



f 

U -ex 



da i 

L7Z 'a=aQSG 



(14) 



As expected, one obtains numerically {Iqsg = 1.0 = 2y6pM- 
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